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bJO ABSTRACT 

We use dimensional regularization to compute the one loop quantum gravi- 
tational contribution to the vacuum polarization on flat space background. 
Adding the appropriate BPHZ counterterm gives a fully renormalized result 
which we employ to quantum correct Maxwell's equations. These equations 
are solved to show that dynamical photons are unchanged, provided the free 
state wave functional is appropriately corrected. The response to the instan- 
taneous appearance of a point dipole reveals a perturbative version of the 
long- conjectured, "smearing of the light-cone". There is no change in the 
far radiation field produced by an alternating dipole. However, the correc- 
tion to the static electric field of a point charge shows strengthening at short 
distances, in contrast to expectations based on the renormalization group. 
We check for gauge dependence by working out the vacuum polarization in 
a general 3-parameter family of covariant gauges. 
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1 Introduction 

Electromagnet ism provided the first example of a relativistic, unified gauge 
theory. Later on, it was quantum electrodynamics (QED) which produced 
the first quantitative successes in the struggle to understand interacting quan- 
tum field theories. It is therefore natural to wonder what electromagnetism 
can tell us about quantum gravity. 

Efforts along these lines date back more than half a century, and were at 
first concerned with a phenomenon termed, "smearing of the light-cone" [I]. 
The idea is that quantum gravitational effects might soften the divergences 
of other quantum field theories because those divergences are associated with 
the singularities all propagators develop for null separations [2J, 



Here iA[g](x; x') is the scalar propagator in the presence of a general metric 
background g^ u , and a[g](x;x') is ~ times the square of the geodesic length 
from x M to x" 1 in that metric. Although the propagator is a well-defined 
distribution — its 4-dimensional integral against a test function converges 
- powers of it are not. That is is why there is a quadratic ultraviolet 
divergence in the two loop "setting sun" contribution to the A0 4 self-mass- 
squared depicted in Fig. [TJ 



For fixed the singularity occurs at different points x ,fM as the metric g^ v 
is varied. Quantizing gravity entails functionally averaging (J2J over metrics, 
and this might be expected to reduce or eliminate the singularity. 



Figure 1: A two loop contribution to the self-mass-squared in A0 4 theory. 
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Of course smearing the light-cone would distort the propagation of light, 
and even a tiny angular deviation might show up over cosmological distances. 
Metric fluctuations can also induce luminosity and redshift variations, as well 
as spectral line broadening and angular blurring. These effects have typically 
been studied indirectly, by working out the scattering of light by a classical 
metric [3J. A significant example is the computation Sachs and Wolfe made 
of the anisotropies cosmic microwave photons acquire in propagating through 
primordial cosmological perturbations [I] . More recent treatments have been 
given of the effects of various ensembles of gravitons [5]. (It is but a short 
step to adjust the distribution to reproduce quantum 0-point fluctuations, 
but this does not seem to have been done.) And there has been much recent 
interest in metric fluctuations engendered by quantum oscillations of matter 
fields [6], a representative diagram for which is given in Fig. [2j 

Quantum gravity also affects electrodynamic forces. Radkowski seems 
to have performed the first computation of the one loop correction to the 
Coulomb potential of a point charge [TJ. A somewhat different result was 
inferred from the scattering of charged scalars [8J. Both of those studies found 
that quantum gravity strengthens the electrostatic force at short distances. 
Much recent interest attended the proposal of Robinson and Wilczek that 
the effect goes the other way, based on renormalization group flows [9]. 

A widely perceived impediment to all of these studies is the fact that 
general relativity plus electromagnetism is not perturbatively renormalizable, 
even at one loop order [10J. One consequence is that we cannot compute ev- 
erything reliably, but that does not mean we are unable to calculate anything. 
One must simply view quantum gravity as a low energy effective field theory 
whose divergences are absorbed using BPHZ (Bogoliubov, Parsiuk, Hepp and 
Zimmermann) counterterms [UJ. In this context, loops of massless gravitons 
and photons engender nonlocal and ultraviolet finite contributions to the 
effective action which are unique predictions of the theory that cannot be 



Figure 2: Lowest order gravitational effect on light due to quantum fluc- 
tuations of matter. Photon lines are wavy, graviton lines are winding, and 
matter lines are solid. 
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changed by its ultraviolet completion [12J. This is why Bloch and Nordsieck 
were able to resolve the infrared problem of QED [T3], long before that the- 
ory's renormalizablility was understood. It is also why Weinberg was able to 
derive a similar resolution for the infrared problem of quantum gravity [14J, 
and why Feinberg and Sucher were able to use Fermi theory to compute the 
long range force engendered by the exchange of massless neutrinos |15j . 

One is of course free to criticize the study of quantum gravity plus elec- 
tromagnetism on the grounds that the predicted effects are too small to be 
observable. Our own interest in the subject derives from its relevance as the 
flat space correspondence limit of the regime of primordial inflation, during 
which quantum gravitational effects are not unobservably small. Indeed, the 
gravitational response to quantum fluctuations of matter [T6] has been re- 
solved [17], and the corresponding fluctuations of gravitational radiation [18] 
may soon be detected [19|. 

Whether during primordial inflation or on flat space background, the 
proper vehicle for studying quantum distortions of electrodynamics is the 
quantum-corrected Maxwell equation. One gets this by first computing 
the "vacuum polarization" zpTF^x; x'), which is the one-particle-irreducible 
(1PI) 2-point function for the photon. This is then used to quantum-correct 
Maxwell's equation, 



This framework has been employed to infer the effects of inflationary charged 
scalar production on photons [20J, and on electrodynamic forces [21]. The 
purpose of this paper is to facilitate a similar study of the effects of inflation- 
ary graviton production by first working out the flat space correspondence 
limit. An example of the utility of this exercise is the recent examination of 
the effects of inflationary scalars on gravitons [22], for which the flat space 
limit [22] provided crucial guidance in dealing with the vastly more compli- 
cated graviton self-energy [24] that pertains during primordial inflation. 

The quantum gravitational contribution to the one loop vacuum polariza- 
tion is derived in section 2. Section 3 solves the quantum-corrected Maxwell 
equation (J3J) for photons, for the instantaneous creation of a point dipole, for 
an alternating point dipole, and for a static point charge. The issue of gauge 
dependence is discussed in section 4, and our conclusions comprise section 5. 




(3) 
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2 One Loop Vacuum Polarization 



The purpose of this section is to compute the renormalized, one loop con- 
tribution to the vacuum polarization from quantum gravity on flat space 
background. We begin by presenting the necessary Feynman rules. Then 
we use them to compute the dimensionally regulated result. By a process 
of successive partial integrations this is expressed as a divergent, local term 
- which is canceled by a BPHZ counterterm — plus the finite, nonlocal 
contribution which constitutes the renormalized result. 



2.1 Feynman Rules 

Our total Lagrangian contains three parts, 

£ — £gr + £em + £bphz ' (4) 

These are, respectively, the Lagrangians of general relativity, electromag- 
netism and the BPHZ counterterm required for this computation, 

C GR = ^RV^, (5) 

£em = —F^F^g^y^, (6) 

Abphz = C 4 D a F^DpF pa g a PgWg™^g- . (7) 

We employ a D-dimensional, spacelike metric g pu , with inverse g^ v and de- 
terminant g = det(g pu . Our affine connection and Riemann tensor are, 

1 

T 

= Q vp -Q vp +vp r a -VP r a (a) 



r V = o9 Pa \duga^ + d!,gua-dagnu\ , (8) 



Our Ricci tensor is R pu = R p ^ pu and the associated Ricci scalar is R = 
g^R^y. The electromagnetic field strength tensor and its first covariant 
derivative are, 

F pv = d^A v -d u A p , (10) 
F> F = f) p _ft f — T 7 F fill 

We define the graviton field h pv { the difference between the full 
metric and its Minkowski background value r]^, 

g^( x ) = Vnv + kV( z ) > ( 12 ) 
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where k 2 = 16nG is the loop counting parameter of quantum gravity. We fol- 
low the usual conventions whereby a comma denotes ordinary differentiation, 
the trace of the graviton field is h = rj^h^, and graviton indices are raised 
and lowered using the Minkowski metric, h^ u = rf p h pv and h pv = rf p r\ va h pa . 
After extracting a surface term the gravitational Lagrangian can be written 
as, 



£gr — Surface 



x< ^hctp^hyaji hap^phu^^ -\- h a p^ph pL , tU ^hctp^hpcr^ / . (13) 



The quadratic part of the invariant Lagrangian is, 



C 



(2) 
GR 



2 ~ ' m°".p 2 
We fix the gauge by adding, 



£ G Rfix = -^F^ 



^h^h^-lh^h^h^h^-h^hp^ . (14) 



^ = ^(v*-^w)- ( 15 ) 



The resulting graviton propagator can be expressed in terms of the massless 
scalar propagator iA(x; x'), 



fivApu 



2 

(x;x') = ^ri^ri^-jj-^ri^rip^iAfax') . 



(16) 



The spacetime dependence of the scalar propagator derives from the Lorentz 
interval Ax 2 (x;x'), 



Ax 2 (x; x') 



x—x 



-{it-t'l-ie)' 



r(——i) 1 5- 

iM.r:.r>) 7 



47TT ^AX 2 



The quadratic part of the electromagnetic action is, 

C EM = ~dpA u d»A» + ^(dpAn 2 ■ 
We fix the gauge by adding, 



(17) 
(18) 



C EMf i x = --(dpAn 2 



(19) 
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Figure 3: Graviton contributions to the one loop vacuum polarization. 



The associated photon propagator is, 



(x;x') = rj IJ ,piA(x;x') . 



(20) 



Electromagnetic interaction vertices descend from the second variational 
derivative of the action, 



5 2 S 



EM 



g(x) g KX (x)g^(x)-g Kp (x)g^(x) d x 6 D (x-x' 



(21) 



5A^x)5A p (x>) 

The necessary vertex functions are obtained by expanding the metric factors, 

^/~~9 {q kX 9^ p ~ Q Kp = rj KXr q pp — rj Kp rj XfJ ' 

+kV ppkX ^K p + K 2 U ppnXa ^ & h a ph 1& + 0(k 3 ) . (22) 

The 3-point and 4-point vertices are, 

yHpKXa/3 _ ^a/3^K[A^p]^_|_^a)[^ft][p^A](/3 ^3) 
_ 1 1 

JJ-fJ,pK\ctfl~f8 



^[-V^pIm _|_ k][p X](d 



_|_^7<5^a) IPyj*] [P^ A ] iP _|_ ■n^-nP) [ A ^P] (7„*)A* _|_^ k (7^5) [A p] {a P)n _|_ n(a ^){-y^S)[\ p]p, 
_|_^ K (7^<5)(a^/3)[A^p]p _|_ ^K[A^p](a^)(7^(5)At _|_ ^K[A^p](7^<5)(a^/3)p ^ ^4) 

Note that parenthesized indices are symmetrized, whereas indices enclosed 
in square brackets are antisymmetrized. 



2.2 Dimensionally Regulated Result 

The three one loop diagrams which contribute to the vacuum polarization 
are depicted in Fig. [3J They can each be expressed using the notation of the 



previous section. The left hand diagram is, 



"IT 



3pt 



[x: x = w 



) 2 d K d'Jv^ KXa ^ 



(x-x')V" a ^ 5 d x d'J 



P A a 



[x: x 



(25) 

Substituting expressions (ITS|) . (12"0)) and (I2"3"j) . acting the inner derivatives and 
performing the inner contractions gives, 



3pt 



167T 



^D[Ax» Ax p r) V(T -Ax p Ax a rj pv ] D[Ax^ Ax v -q^ - Ax p Ax v ri^U 



Ax 



2D 



Ax 



2D 



The next step is to act the outer derivatives, at which point we can extract 
a manifestly transverse form, 



3pt 



2 tl fn ,vn o^nJ P+l)r 2D^Ax" 



/t 2 r 2 (f-i) (d-3)(d-2) 2 d 

" 16V 5 2(D-1) 
The middle diagram of Fig. 1 is, 



rf v d 2 -d il d h 



Ax 2D + 2 
1 



4pt 



(x;a/) = m 2 <9 K jf/^ KAa/37<5 z 



Ax 



(x; x)d\5 D (x — x' 



2D-2 



(27) 
(28) 



(29) 



This diagram vanishes because the coincidence limit of the massless scalar 
propagator in flat space is zero in dimensional regularization, iA(x;x) = 0. 
The diagram on the right of Fig. 1 is, 



p W 



ctm 



[x: x 



') = i4Cjrj pu d 2 -d p d u )d 2 5 D (x-x' 



(30) 



2.3 Renormalization 

To renormalize ( 128|) we must first localize the ultraviolet divergence so that 
it can be subtracted by the counterterm (I3"0"j) . This process of localization 
is accomplished by first partially integrating the factor of l/Ax 2D_2 in (12 8p 
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until the remainder is integrable [21] • In dimensional regularization the steps 



arc 



1 



cf 2 



1 



Ax 2D - 2 2(D-2) 2 Ax 2D ~ 4 4(D-2) 2 (D-3)(D-4) Ax 2D ~ 6 
Next we add zero in the form 



d 2 



i 



Ax D ~ 2 r(f-i) 



5 D (x-x') 



(31) 



(32) 



Adding 

d 2 



to the key part of (13"Tj) in a dimensionally consistent way gives, 



D-4 Ax 



2D- 



iATTT /J, D - 4 5 D (X-X') 



2 



D-4 



r 



D 



D~4 



D-4\ Ax 2D ~ 



Ax D ~ 



i4%% fi D - 4 6 D (x-x') <9 2 / m(/i 2 A 
D-4 ~2~\ A^ 2 



r(f-i; 



0(D-4) 



(33) 
(34) 



Substituting ( 13T|) and ( l34"l) into ( 128]) results in the desired localized divergence, 

D n D 



3pt 



x; x 



z/t 2 r(f-r 



+ 



47T 2 
hi 2 



8(D-l)(D-4) 



tiV d 2 -& 1 d v d 2 S D (x-x' 



192n 4 



fMV d 2 -d fl d u }d /L 



Jln(/x 2 A; 



+ 0{D-4) . (35) 



The local divergence of expression fl35l) will be completely canceled by 
the counterterm fl30l) if we make the choice, 



C 4 



« 2 r(f-i) 



16ttt 8(D-l)(D-4) 



(36) 



We can then take the unregulated limit (D — > 4) to obtain the fully renor- 
malized graviton contribution to the one loop vacuum polarization, 



(x; x') 



IK 



192tH 



^Q 2 _ «9 4 / M^ 2 ^ 2 ) 
-I | Ax 2 



(37) 



Note that the ambiguity regarding the finite part of the counterterm is re- 
flected in the dimensional regularization scale /z. 
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3 Quantum Corrected Maxwell Equations 



The purpose of this section is to use our one result (1371) for the one loop 
vacuum polarization to quantum correct Maxwell's equations, and then in- 
fer quantum gravitational corrections to electrodynamics by solving these 
equations. We begin by deriving the causal effective field equations of the 
Schwinger-Keldysh formalism. Subsequent subsections solve these equations 
perturbatively for the special cases of free photons, a point dipole pulse, an 
alternating point dipole, and a static point charge. 

3.1 Schwinger-Keldysh Formalism 

We come now to the question of what to use for the vacuum polarization 
[TP] (a?; x') in the quantum corrected Maxwell equation (J3J). It might seem 
obvious that the in-out result (l3Tj) we have just derived should be used, but 
that would lead to two problems: 

• Causality — The in-out vacuum polarization ( 137|) is nonzero for points 

which lie in the future of or at spacelike separation from it; and 

• Reality — The in-out vacuum polarization ( l3Tj) is not real. 

One can get the right result for a static potential by simply ignoring the 
imaginary part [8] , but circumventing the limitations of the in-out formal- 
ism becomes more and more difficult as time dependent sources and higher 
order corrections are included, and these techniques break down entirely for 
the case of cosmology in which there may not even be asymptotic vacua. 
Note that there is nothing wrong with the in-out vacuum polarization (I37j) ; 
it is exactly the right thing to correct the photon propagator for asymptotic 
scattering computations in flat space. The point is rather that employing 
fl37|) in equation ([3]) fails to provide a set of field equations with the same 
scope and power as the classical Maxwell's equations. 

The more appropriate field equations are those of the Schwinger-Keldysh 
formalism. This technique provides a way of computing true expectation 
values that is almost as simple as the Feynman diagrams which produce in- 
out matrix elements [27]. We shall develop the Schwinger-Keldysh rules in 
the context of a scalar field ip(x) whose Lagrangian (the space integral of its 
Lagrangian density) at time t is L[ip(t)]. Suppose we are given a Heisenberg 
state whose wave functional in terms of the operator eigenkets at time to 
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is \&[<^(io)], and we wish to take the expectation value, in the presence of this 
state, of a product of two functionals of the field operator: A[ip], which is 
anti-time-ordered, and B[(p], which is time-ordered. The Schwinger-Keldysh 
functional integral for this is [28], 



The time t\ > to is arbitrary as long as no operator in either A [ip] or B[ip] is 
evaluated at a later time. 

The Schwinger-Keldysh rules can be read off from its functional represen- 
tation fl38|) . Because the same field operator is represented by two different 
dummy functional variables, (p±(x), the endpoints of lines carry a ± polar- 
ity. External lines associated with the anti-time-ordered operator A [cp] have 
the — polarity whereas those associated with the time-ordered operator B[ip] 
have the + polarity. Interaction vertices are either all + or all — . Vertices 
with + polarity are the same as in the usual Feynman rules whereas ver- 
tices with the — polarity have an additional minus sign. If the state \^f) is 
something other than free vacuum then it contributes additional interaction 
vertices on the initial value surface [29] . 

Propagators can be ++, H — , — h, or . All four polarity variations 

can be read off from the fundamental relation (138]) when the free Lagrangian 
is substituted for the full one. It is useful to denote canonical expectation 
values in the free theory with a subscript 0. With this convention we see that 
the ++ propagator is just the ordinary Feynman propagator, 



iA ++ (x; x' 



n 



T(tp(x)(p(xf] 



n 



iA(x; x') 



(39) 



where T stands for time-ordering and T denotes anti-time-ordering. The 
other polarity variations are simple to read off and to relate to the Feynman 
propagator, 



iA_ + (x; x 



iA + _(x; x') = (tt 
iA__(x;x')= (ft 



ip(x)ip(x' 
<f(x')<f(x 
T(<f(x)(p(x' 



n)=6(t-t')iA(x;x')+6(t'-t)[iA(x;x')\ , (40) 
tt)= 9{t-t') [iA(x] x')] *+6(t'-t)iA(x; x r ), (41) 



n 



iA(x: x') 



(42) 
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In our case, both the photon and the graviton propagators depend upon the 
massless scalar propagator (1171) . which is a function of the Lorentz interval 
Ax 2 (x;x'). It follows from relations f l40H42|) that the various Schwinger- 
Keldysh propagators can be obtained by making simple replacements for the 
Lorentz interval, 





(x; 


X') EE 


x — 


x 1 


Ax 2 + _ 


(x; 


X') EE 


x — 


x' 


Ax 2 + 


0; 


X') EE 


x — 


x' 


Ax 2 __ 


(x; 


X') EE 


x — 


x' 



2 - c 2 (\t-t'\-ie) 2 



2 - c 2 (t-t'+ie] 2 



l {\t-t'\ + 



if. 



(43) 
(44) 
(45) 
(46) 



Because each external line can be either + or — in the Schwinger-Keldysh 
formalism, every 1PI N-point function of the in-out formalism gives rise to 
2 N 1PI N-point functions in the Schwinger-Keldysh formalism. For every 
classical field <p(x) of an in-out effective action, the corresponding Schwinger- 
Keldysh effective action must depend upon two fields — call them 4>+(x) and 
<p-{x) — in order to access the appropriate 1PI function [30J. For the scalar 
paradigm we have been considering the 1PI 2-point function is the scalar 
self-mass-squared, M 2 ± (x;x'), and the effective action takes the form, 



0_] = S[<f> + ] - S[4>_] 



- \d\ ld 4 x' 



x 



<p+(x)Ml + (x; x')(j) + (x') + (p+(x)Ml_(x; x')(f)_(x') 
+(p_(x)M 2 + (x; x')<p + (x') + (j)_(x)M 2 _(x; x')cj)_{xf) 



0(4), (47) 



where S is the classical action. The effective field equations are obtained by 
varying with respect to + and then setting both fields equal [30] . 



5T\ 



<t>±-- 



O 2 



■m 



x)- d 4 x' Ml + (x;x')+Ml_(x;x') 4>{x') + 0{4> 



(48) 

The two 1PI 2-point functions we would need to quantum correct the lin- 
earized scalar field equation are M 2 + (x; x') and M 2 _(x; x'). Their sum in (jl8"]l 
gives effective field equations which are causal in the sense that the two 1PI 
functions cancel unless x' M lies on or within the past light-cone of x M . Their 
sum is also real, which neither 1PI function is separately. 
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From the preceding discussion it is apparent that we wish to make the 
following substitution in equation (jSJ), 



(x; x') 



(x: x') + 



(x; x') 



(49) 



where we can read off the appropriate Schwinger-Keldysh vacuum polariza- 
tion from expression (131 



±- L1 ± 



(x; x') 



(±)(±)« s 



1927T 4 

±)(±)m 2 



" u d 2 -d fl d l 







I ln(/i 2 A 



(50) 



^ 2 -a^]a 6 {ln 2 (/i 2 Ax 2 t± )-21n(/i 2 A4 ± )} . (51) 



1536tt 4 

Now define the temporal and spatial intervals as, 

At = t-t' , Ar=||f-f||. (52) 
that differences of logarithms of the 



It is apparent from expressions 
the ++ and H — intervals give, 



ln(n 2 Ax 2 ++ ) - \n(n 2 Ax 2 + _) = 2iri6(At-Ar 



ln 2 (/i 2 Aa;; + ) -\n A (^Ax 2 + _) = Am6 (At-Ar In //(At 2 -Ar 



(53) 
(54) 



Hence the vacuum polarization which belongs in equation ([3]) is, 

.2 



(x; x) 



K 



384vr 



x 



r] »v d 2_ d H d v 

a 6 |^(At-A 



In 



fi 2 At 2 -Ar 2 



+ 0(k 4 ) . (55) 



3.2 Photons 

Expression (155]) gives all the one loop contributions which derive exclusively 
from interaction vertices, but there are also contributions from perturbative 
corrections to the initial state wave functionals. In the scalar functional 
integral ( 138]) these wave functionals are \I/[y2 + (to)] and \l/*[<£_(£o)]; for gravity 
plus electromagnetism they would be functionals of and h^, evaluated 
at the initial time. 
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Each state wave functional can be expressed as the wave functional of 
free vacuum times a series of perturtbative corrections, 



V[A, h] = V [A, h] x \l + o(KhA 2 ) 



(56) 



It is straightforward to show that the free vacuum contribution is what fixes 
the real part of the propagator in the functional formalism [31] . If there were 
no perturbative state corrections then merely employing the correct propa- 
gators would completely account for the state wave functionals. However, 
there must be perturbative state corrections because free vacuum cannot be 
the true vacuum state of an interacting quantum field theory. 

Perturbative state corrections manifest as new interactions on the initial 
value surface [28]. When the initial value surface is in the asymptotic past (or 
the asymptotic past and future for in-out matrix elements) these interactions 
have no effect on operators at finite times. However, they can be important 
when the initial value surface is at a finite time, as it must be in cosmology. 
The first correction relevant for a massless, minimally coupled A0 4 theory has 
recently been worked out on de Sitter background [29J. In this case the initial 
state correction is necessary to make the linearized effective field equation 
well defined at the initial time (32], and to eliminate an infinite series of 
rapidly redshifting terms from the two loop expectation value of the stress 
tensor [26] . 

We shall assume that the missing state corrections exactly cancel the 
surface terms which arise when (1551) is partially integrated. To see what this 
entails, first note that all orders of the "pure-vertex" part of the vacuum 
polarization take the manifestly transverse form, 



The partial integration we have in mind concerns the quantum correction to 
Maxwell's equation, 



In the Schwinger-Keldysh formalism the ++ and H — contributions exactly 
cancel on the future temporal surface, as well as on the surface at spatial 



U u ](x;x') = (7]^d 2 -d^d u )u(x-x') . 



(57) 




d 4 x'U{x-x')d' v F VfM (x') + Surface Terms . (58) 
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infinity. Hence the only surface terms come from the initial time. Of course 
this is also true of perturbative state corrections. We assume that the two 
contributions exactly cancel, so that the full, quantum-corrected Maxwell 
equation is, 

d v F vlx {x) + Jd 4 x'U(x-x') d ,u F^(x') = J^x) . (59) 

We are finally ready to consider the case of free photons, which corre- 
sponds to J^(x) = 0. Note from equation (1511 that these obey d v F y ^{x) = 0, 
the same as in the classical theory. One might worry about the potential 
for solutions of the form d u F u/1 (x) = S^x), where S^x) obeys the integral 
equation, 

S„(x) + Jd A x'Xl{x-x')S^x') = . (60) 

However, an effective field equation such as fl59"|) can only be used to pertur- 
batively correct classical solutions [33], which means we must exclude any 
such solutions. Hence we conclude that quantum gravity on flat space back- 
ground makes no correction to free photons at any order, except for possible 
field strength renormalization. 

3.3 Instantaneously Creating A Point Dipole 

The charge density of a static point electric dipole p at the origin is p = 
—p ■ V5 3 (x). We might imagine creating such a dipole at the instant t = 
by separating the charges in a very small, neutral particle such as a neutron. 
The conserved 4-current associated with this event is, 

J°(t,x) = -6{t)p-V5 3 {x) , J\t,x) =p i 5{t)5 3 {x) . (61) 

The response of the magnetic field provides a good perturbative illustration 
of the smearing of the light-cone which was conjectured so long ago pp. 

Before proceeding it is desirable to reorganize equation (I59"|) in two ways. 
The first has to do with the limitation inherent in only possessing the first 
order term in the loop expansion of U(x — x'), 

U(x-x') = nW(i-i') + U^\x-x') + 0{k 6 ) . (62) 

Of course this means we can only infer the one loop correction to the field 
strength, so we may as well expand it, 

F, v {x) = FjS(x) + F${x) + Fff(x) + 0(n e ) . (63) 
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Substituting (1521) and (1551) in the quantum-corrected Maxwell equation (1551) 
and segregating different orders of k 2 produces the hierarchy, 



Jy.{x) 



d 4 x'U(x-x')J^x') 



(64) 
(65) 



and so on. Note the classical source J M (x) is Oth order. 

The second reorganization concerns deriving the field strength directly, 
without constructing the vector potential. Consider taking the curl of the 
classical Maxwell equation, 



patiuQ gap 



d 2 6 p °» v d lx A v = e^d^Jv 



d 2 F, 



8 T — 8 J 



Combining this with ( I641l65|) implies, 

& p(o) _ f) j _ f) j 

d 2 F^ = d^-hJP. 



(66) 

(67) 
(68) 



Now recall that our one loop current density can be expressed as the d'Alem- 
bertian of something, 



WW 



Jd l x'n ( > 1 \x-x')J ll (x') 



d x 



, fln(fi 2 Ax 2 ++ ) \n(n 2 Ax 2 + _) 
{ Ax 2 ++ Ax 2 + _ 

dV 9 (At - Ar) j In [/i 2 (At 2 - Ar 



IK 2 8^ 
1927T 4 
K 2 8 6 

384tt 3 

Comparison of (168]) with ( 170|) or (}7T|) implies a result for the one loop field 
strength, up to possible homogeneous terms, 



(69) 

W) , (70) 
-lU(x'). (71) 



, j ln(/! 2 Ax 2 + ) ln(/i 2 Aa 
j Ar 



2 

++ 



Ax 2 



>J u ](x') 



(72) 



28[Jd 4 x 



28 [ Jd 4 x'6(At-Ar)hn[ii 2 (At 2 -Ar 2 )}-l\j iy] (x') . (73) 



in 2 8 2 

1927T 4 

~ 384tt 3 

We are now ready to specialize to the current density ( IBTj) of an instanta- 
neously created dipole. Substituting in (jMj) and specializing to purely spatial 
indices gives, 



&F§\x) = (d iPj -d j p t )6 4 
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(74) 



The solution can be expressed in a convenient form by noting the D = 4 
dimensional version of relation 



Wt^1=4^ 4 (*-*') , W-i^)=0. (75) 



Hence we have, 



where x /At = is understood. Now write out the two intervals, 

Ax 2 ++ = r 2 - t 2 + e 2 + 2e|t|i , (77) 
Ax 2 + _ = r 2 - t 2 + e 2 - 2etz . (78) 

Combining these relations with the Dirac identity results in the familiar form 
for the Lienard-Wiechert potential, 

F^\x) = -^(d iPj -d jPi )e(t)5(r 2 -t 2 ) . (79) 

The most convenient form for the one loop correction is f!72|) , 

(rln 8^ ^- P n (/^+ + ) ln(/i 2 Ax 2 J | 

= (^-^)1^{a^;- a^} ' (81) 

-(^)(^){^}- 

Adding the one loop magnetic field (!82|) to the tree one (179|) leads to an 
interesting form, 

= -^(^-a,p^(t)5(r 2 -t 2 - I ^)+C»(K 4 ). (84) 

It would therefore be fair to say that, by time, t the signal has reached a 
distance r slightly outside the classical light-cone, 

r 2 = t 2 + ^ + CV) . (85) 
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Although intriguing, the super-luminality we have just found is unob- 
servably small. In particular, it cannot serve as any sort of explanation for 
the OPERA result [34J. It also isn't cumulative, so looking at cosmological 
sources makes the effect no larger. Another thing our effect fails to do is break 
Lorentz invariance, as could have been predicted from the fact that pertur- 
bative quantum gravity provides no mechanism for spontaneously breaking 
this symmetry. Instead of signals propagating along the classical light-cone 
fi^x^x" = 0, they now propagate along r] IMU x fJ, x 1 ' = So it is not that that 
the speed of light or the dispersion relation has been changed. 

Of course theories with a nonlinear kinetic operator can show super- 
luminal propagation even classically [35J. Our effect is different in that it 
arises from quantum fluctuations of the metric operator which sets the light- 
cone. One interpretation for the net super-luminal propagation is that there 
is more volume outside the classical light-cone than inside. This might be 
checked by extending the graviton expansion of the volume of the past light- 
cone one order higher than in [36] and then computing its expectation value. 
If the one loop correction is positive then our conjecture is verified. Note 
also that this check would be independent of the choice of gauge because the 
volume of the past light-cone is a gauge invariant operator. 

There have been many claims of super-luminal propagation from quan- 
tum electrodynamics in nontrivial geometries [37, 38J. Our result is different 
in that it occurs in flat space, and is due to fluctuations of the metric op- 
erator, rather than of some matter field. We also doubt that the earlier 
claims result from true super-luminal propagation. One cannot compute the 
vacuum polarization produced by fermions in an arbitrary geometry because 
the fermion propagator is not known for general metric. What was done 
instead is a derivative expansion. This should be valid for low energy ef- 
fective field theory; for example, it should give correct results for the phase 
velocity of some continuous, low frequency signal. However, demonstrating 
true super-luminality requires following the propagation of a pulse, and the 
high frequency modes which are essential for this are not correctly treated by 
derivative expansions. In fact the Schwinger-Keldysh formalism [27] implies 
there cannot be super-luminal propagation from the fermionic contribution 
to vacuum polarization. 
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3.4 An Alternating Point Dipole 

The 4-current associated with an alternating point dipole is, 

J°(t,x) = -p-V5\x)e~ iut , J l (t,x) = -iup l 5\x)e- wt 



(86) 



To find the quantum correction to the current we employ the same expansion 
technique used in the previous section where the first order correction is 
defined as, 

Jji)(z) = --^jd 4 x'6(At - Ax) {in [/i 2 (At 2 -Ax 2 )] -ljj'V) • (87) 

We can evaluate this integral by rewriting x and the differential operators as 
x = \{x + 1) + \{x — t) and d 2 = \{d x — d t ){d x + d t )x. Thus we come to the 
convienient form, 

<9 4 = ^(d x -d t )\d x +d t )\x+t)+^(d x -d t )\d x +d t )\x-t) . (88) 



By subsituting (|88|) for d 4 in (157)) and applying the zeroth order currents (1861) 
we find the one loop currents to be, 



jL{t,x) = d 2 



Gv-V v icu 1 
+ 



[ 67T 2 

J iup l G 
| 6n 2 



x 2 x 3 



icu 1 

X 2 X 3 



-iu>(t—x) 



-iu){t—x) 



(89) 
(90) 



From flBTj) we see that the zeroth order field strengths for this source obey, 



d 2 F$\t,x) 

}2 p(0) 



d F-j } (t,x = —iu diPj — djp. 



5 3 (x)e- iuJt 
S 3 (x)e- iuit 



Applying the Lienard-Wiechert Green's function we find, 



F^(t,x) 



4tt 



[uj 2 pi-dip-V 



IUJ r. 



X 

-iw(t—x) 



X 



(91) 
(92) 



(93) 
(94) 
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From (I68p we see that one loop field strengths follow by simply deleting the 
d 2 from f)89H90p and acting some derivatives, 



iuG 



6?r 2 
iuG 



6n 



u 2 pi-dip-V 



1 + - 



uox 



1+- 



ux 



—iu>(t—x) 

X A 
^—iui(t—x) 

X 2 



Adding the loop correction to the tree results gives, 



F 0i (t,x)= u 2 pi-dip-V 



Fij (t,x)= iu diPj - djPi 



x 



3 — iw(t— x) 

47rx 
Atcx 



X 1 



2iuG 



x 1 



3nx 
2iuG 
3nx 



1+- 



USX 



1 + — 

uix- 



+ 0(G 2 



(95) 
(96) 

(97) 
(98) 



Of course the obvious conclusion is that the one loop corrections have no 
effect in the far field regime, and the near field regime is unobservably close 
to the source. 



3.5 A Static Point Charge 

The charge density of a static point charge q at the origin is p = q5 3 (x). The 
conserved 4-current associated with this source is, 

J"(x) = q6 3 {x)5% . (99) 

Because the \x = component differs from the alternating dipole of the 
previous subsection only by setting uj = and replacing — p-V with q, we 
can read off the one loop current density from (|89[) . 

4 } (t,x) = -^L. (ioo) 

Of course the vector components vanish so we find the correction to the 
Coulomb potential is, 

*< r >-iM i+ !! +o < G2 >)- (101) 

Our result ( 1101 p agrees with that found in 1970 by Radkowski j7]. The one 
loop correction that Bjerrum-Bohr inferred from the scattering of charged, 
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Figure 4: Vertex correction included (along with many other diagrams) in 
the Bjerrum-Bohr result [8], but not in either our result fllOip or that of 
Radkowski [TJ. Charged scalar lines are solid with an arrow, photon lines are 
wavy and graviton lines are winding. 

gravitating scalars differs from what we got by a factor of nine [8]. Part of 
this discrepancy may be due to different sources; Bjerrum-Bohr considered a 
charged scalar whereas we used a point particle with worldline x M ( r )) 

Vint = -rnyJ-g^( X {t))x*{T)x v {T) + gX M (r)A M ( X (r)) . (102) 

However, we believe the largest part of the discrepancy arises from Bjerrum- 
Bohr having implicitly included corrections to the current density like the 
diagram depicted in Fig. HI We could have and should have done this, but 
we will see in the next section that it would only have altered all our one 
loop field strengths by an overall constant. 

Radkowski [7J, Bjerrum-Bohr [8] and we all agree that quantum gravity 
strengthens the electromagnetic force at one loop. The opposite conclusion 
seems to arise from computations of the quantum gravitational contribution 
to the electromagnetic beta function [91 [39j HQ]. These show that quantum 
gravity decreases the electromagnetic coupling constant at high energy scales. 
That would normally be assumed to mean that quantum gravity weakens 
the electromagnetic force at short distances, but it is well to keep in mind 
that the beta function is not directly observable. The observable thing is 
the strength scattering between charged particles, and the Bjerrum-Bohr 
computation shows that one loop quantum gravity effects weaken this, rather 
than strengthening it. 
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4 Gauge Dependence 



The purpose of this section is to examine how the results of the previous 
section depend upon our choices of the gravitational gauge fixing term ([15!) 
and the electromagnetic gauge fixing term ( Tl9l) . We begin with some general 
considerations which reduce the issue to a single proportionality constant. 
The graviton and photon propagators are then worked out for a general 3- 
parameter family of covariant gauges. Although one of these parameters 
drops out, the other two can change the proportionality constant all the way 
from minus infinity to plus infinity. We close by exploiting the gauge inde- 
pendent result of Bjerrum-Bohr to argue that this seeming gauge dependence 
may cancel out if quantum gravitational corrections to the current density 
are included. 



4.1 General Considerations 

Note from expression (|2^|) that the vacuum polarization is transverse on each 
of its two indices as a trivial consequence of the antisymmetry of the vertex 
function on its first and third indices, 

ytipK\a/3 _ __yKpii\a/3 (103) 

This is completely without regard to the gauges employed to define graviton 
and photon propagators. Suppose we now restrict attention to gauges which 
preserve Poincare invariance. Because the Lagrangians ([MI]) and the back- 
ground are also Poincare invariant, the vacuum polarization must inherit this 
symmetry. Then dimensional analysis, transversality and the standard k 2 of 
a one loop quantum gravity result, together imply a form like that of ( 128]) . 



iJ -3pt 



(x-x) = -Constant x k 2 ^ v d 2 -d l "d v A on - . (104) 

L J /\x 2U 1 



However, the constant prefactor might be gauge dependent, and that same 
gauge dependent constant would multiply all of our one loop corrections. 

It is useful to begin at a somewhat earlier point. If different - - but 
Poincare invariant — graviton and photon propagators had been employed 
in expression (|2"5|) . then the combination of Poincare invariance, dimensional 
analysis and the algebraic symmetries of the vertex function and the prop- 
agators imply that the result can be expressed in terms of two constants A 
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Cj(D, a, b) 




1 


1 




2 


2 

D-2 


VfJ.i'Vpcr 


3 


4(6-1) 
(D-2)(6-2) 


dpda | 

'/M^ 9 2 1 <g2 '/per 


4 


a — 1 


* d 2 


5 


4o(6-l)(6-3) 8(6-1) (6+-D-3) 




(6-2) 2 (6-2) 2 (D-2) 


a 4 



Table 1: Coefficient functions Ci(D, a, b) and the tensor differential operators 
[ftvTpo] for the graviton propagator fill II) defined with the general gauge fixing 
functional (11 101) . 
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D 



2P-1) 
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(105) 



Ax 2D- 



(106) 
.(107) 



We can therefore identify the proportionality constant in fjl04j) as, 
Constant = E!<izi> x (0-2)[ZM-2(B-l)B] 



16tt d 

r 2 (f-i) i /D—2 



16tt 



d 



x 



2(£>-l) 
2VD-1V 



(108) 
(109) 
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4.2 General Covariant Gauges 



The most general Poincare invariant extension of the graviton gauge fixing 
functional ffT5]) depends upon two parameters a and 6, 



2a V J-*** 



•^GRnew 

The associated propagator is jUJ, 



T = n pa ( h h 



Anew 



J2Ci{D,a,b) x L7/J x tA{x;x' 



1=1 



;no) 



in) 



where the coefficient functions Ci(D, a, b, ) and the tensor differential opera- 
tors [nvT^] are given in Tabled) The propagator can be given a more reveal- 



ing expression using the transverse projection operator = rj^ 



Q2 ) 



Anew 
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^p(p^cr)u- 
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( 6-2 
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<9 2 

<9 4 



.(112) 



Of course the tranverse-traceless term on the first line represents the con- 
tribution from dynamical, spin two gravitons. This term looms large in the 
quantum gravity literature but it is well to recall that it plays no role in the 
solar system tests of general relativity. The phenomenologically more im- 
portant parts of the graviton propagator are those on the second and third 
lines, which mediate the gravitational interaction between sources of stress- 
energy. Note that the longitudinal terms proportional to the gauge parameter 
a would vanish in the exact gauge h v = \h^. 

The most general Poincare invariant extension of the photon gauge fixing 
functional ( TT9|) depends upon a single parameter c, 



£ EMnew 

The associated propagator is, 
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2c 
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[x: x 



T] pa + (c-1) 



d p d a 
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(113) 



;ii4) 
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The longitudinal term proportional to c — 1 can make no contribution to the 
general gauge vacuum polarization f)105|) because the vertex function (I2"5j) is 
antisymmetric under interchange of its second and fourth indices, 

y(j,pK,\a(3 _ _y(j,\KpaP (115) 

It remains to explain how to act the tensor differential operators of Tabled] 
on the scalar propagator ffTTj) . First note that inverse d'Alembertians act on 
l/Ax D ~ 2 to give, 



11 11 



d 2 Ax°- 2 ' 2(D-4)Ax D - 4, 
11 11 



d 4 Ax D - 2 8(D-4)( J D-6) Ax D " 6 ' 
Now just act the derivatives in the numerator to conclude, 



(116) 
(117) 



dpd u . I f {D-2)Ax,Ax u \ . A , 

-jp-iA{x;x) = 2 X }^ n&(x;x) , (118) 

dpd v d p d a 1 J Q(D-2)ri (pv Ax p Ax a) 

+— A M x4 iA(a;; . (120) 



4.3 Gauge Dependent Proportionality Constant 

We are now ready to compute the crucial proportionality constant of rela- 
tion (I104p . Because the gauge dependence of the photon propagator drops 
out, we need only consider the gauge dependence of the graviton propaga- 
tor. Because the graviton propagator (11 111) is a sum of gauge-dependent 
coefficients Ci(D,a,b) times tensor operators [fiuTL], acting on the scalar 
propagator, we may as well work out the result for each tensor operator 
separately. Table [2] presents the coefficients Aj(D) and Bj(D) which were 
defined in relation (I105p . for each of the five tensor differential operators of 
Table [TJ Also given is the contribution of each tensor differential operator to 
the coefficient Cj(D), 

C^D) = D x A T (D) - 2(D-1) x B T (D) . (121) 
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Table 2: Coefficients A/, 5/ and Cj defined in relations (jll)5p and ( I109p . 
under the replacement i[ Ali/ A°° w ](x; x') — )■ [^TL] X iA(x; x') for each of the 
five tensor differential operators defined in Table [D 



C(4, a, 6) = -12a x 1 ^ ^ ; - 7 ^ r — - . (123) 



To recover the full result for C(D, a, 6) defined in relation (I109p we mul- 
tiply each Ci(D) by the appropriate gauge dependent coefficient C(D, a, 6) 
from Table HJ 

5 

C(D, a,b) = Y^ Ci{D, a, b) x d(D) (122) 
i=i 

The formula for arbitrary D is not illuminating, but specializing to D = 4 
gives, 

(36 2 - 126+8) 

(6-2) 2 (6-2) 2 

Our original gauge corresponds to a = b = 1, which gives C(4,a, 1) = 8. 
Hence the various one loop corrections computed in section are valid for 
general a and b if we multiply by the proportionality constant, 

r „ IN C(4,a,6) 3 ( 36 2 - 126 + 8) 1 

It is interesting to note that the gauge independent contribution from dy- 
namical gravitons vanishes in D = 4 dimensions, 

C 1 (D)-C 4 (D) + 2C 5 (D)-^[c 2 (D)-C 3 (D)+C 5 (D) 

_ (D-A){D-2)\D + l){D + 2) 

~ 4(^1) - (i25j 
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It is apparent that the gauge dependent proportionality constant (I124p 
can be made to have either sign by varying the gauge parameter a. Further- 
more, K(a, b) can be made arbitrarily large in magnitude by taking the gauge 
parameter b close to 2. Hence it would seem that our results are completely 
gauge dependent and unphysical. Gauge dependence has also been noted in 
the renormalization group approach [4"2"j |4"5] . 

A moment's thought reveals that all is not lost because the result of 
Bjerrum-Bohr [8] for the quantum gravitational correction to the Coulomb 
potential was derived from the gauge independent S-matrix of scalar QED. 
Roughly speaking, this correction derives from the fact that gravity is sourced 
by the electromagnetic fields of the two charged particles being scattered, and 
this source changes as the particles move with respect to one another. That 
is a real effect, not some gauge artifact. And it is crucially important to 
note that we agree with Bjerrum-Bohr up to a factor of +9. We attributed 
this factor to our having only quantum-corrected the left hand side of the 
operator Maxwell equation, 



The right hand side is also an operator and it must also suffer quantum 
gravitational corrections, one of which is depicted in Fig. HJ By Poincare 
invariance, current conservation and dimensional analysis, those corrections 
must take exactly the same form as we found for the left hand side, up to 
an overall constant. We conjecture that the gauge dependence K(a, b) we 
have just found for corrections to the left hand side is canceled by gauge 
dependence in corrections to the right hand side. If this is correct, then the 
overall, gauge independent correction to the various results derived in section 
3 can be inferred by comparing any one of them with its S-matrix analogue. 
For scalar QED the correction factor would be 9, and it could be computed 
for the point particle source we used. 

The resolution we have just proposed to the gauge issue recalls some old 
work by DeWitt [2] about dependence upon the gauge fixing functionals 
even in the gauge invariant background field effective action V = S + E. 
DeWitt states |45j, "The functional form of £ is not independent of the 
choice of these terms. However, the solutions of the effective field equation 
can be shown to be the same for all choices." At the order we are working 
there is no distinction between our effective field equations and those of the 
gauge invariant background field effective action. (One can see this from the 




(126) 
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transversality of our vacuum polarization.) However, we have just shown that 
De Witt's statement cannot be correct if the source (or the asymptotic field 
strengths for scattering solutions) is not normalized in some physical way. 
Our proposal is that including quantum corrections to the right hand side of 
the equation provides this physical normalization. More work is obviously 
required, in particular an explicit computation of the quantum gravitational 
corrections to the source, but it would be wonderful if solutions to the effec- 
tive field equations could be physically interpreted the same way as classical 
solutions. 

5 Discussion 

We used dimensional regularization to compute the one loop quantum grav- 
itational contribution ( 128|) to the vacuum polarization on flat space back- 
ground. A fully renormalized result (|37|) was obtained by first partially inte- 
grating to localize the ultraviolet divergence and then absorbing it into the 
appropriate BPHZ counterterm (1301) with coefficient (I36p . The Schwinger- 
Keldysh formalism [271 [30] was then employed to reach the manifestly real 
and causal form fl55l) . 

We used f )55|) to solve the quantum corrected Maxwell's equation ([3]) 
for various special cases. Provided the appropriate perturbative corrections 
to the initial state cancel the surface terms involved in reaching the form 
(|59|) . there is no change in the source-free solutions at any order in the loop 
expansion. However, sources induce a variety of interesting effects. 

Probably the most provocative source is the current density (|6Tj) of an 
instantaneously created, point electric dipole. The pulse (ISlj) which results 
in the magnetic field propagates slightly outside the classical light-cone. It 
seems to arise from quantum fluctuations of the metric operator, which are 
isotropic but favor super-luminal propagation because there is more volume 
outside the light-cone than inside. That this sort of thing might occur has 
been realized since the earliest days of quantum gravity [TJ. Our super- 
luminal effect is completely Lorentz invariant, merely changing the charac- 
teristic surface from r] IM ,x^x v = to r\ liV x Pj x v = |p. Despite many claims 
to the contrary, this seems to be the first case of super-luminal propaga- 
tion from a quantum field theory whose classical analogue does not allow 
super-luminal propagation. All previous claims have been based derivative 
expansions [37J 135] , which are perfectly valid for most applications of low 
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energy effective field theory but which incorrectly treat the high frequency 
modes needed to resolve the propagation of a pulse. 

The other interesting source we studied is the response to a static, point 
charge. Our result (llOip for the quantum corrected Coulomb potential agrees 
with what Radkowski found more than four decades ago [7] . It does not agree 
with the Coulomb potential Bjerrum-Bohr inferred from the scattering of 
charged, gravitating scalars [8J, but we all find quantum gravity strengthens 
the electrostatic force at short distances. We believe that the factor of nine 
discrepancy with Bjerrum-Bohr derives from his S-matrix technique implic- 
itly including quantum gravitational corrections to the charge density, like 
the diagram depicted in Fig. |H This should have been done for our point 
particle source, but it would only have changed the one loop field strengths 
by an overall constant. 

Renormalization group analyzes [9~J [39j HQ] seem to provide a more serious 
discrepancy. These studies find that quantum gravity reduces the electrody- 
namic coupling constant at large scales. The usual inference would be that 
quantum gravity weakens the electrostatic force at short distances. However, 
the beta function is not itself observable; several other effects must be com- 
bined to infer the impact of quantum gravity. The S-matrix computation of 
Bjerrum-Bohr should include all of these effects, and it shows that quantum 
gravity strengthens the electrostatic force at short distances. 

Gauge dependence poses a major obstacle to the physical interpretation of 
solutions to the effective field equations. If one restricts to Poincare invariant 
gauge fixing functionals, the only possible change to our vacuum polarization 
f )55|) is rescaling by an overall, gauge- dependent constant. In section 4 we 
considered the most general 2-parameter family of graviton gauges (IllOp , and 
the most general 1-parameter family of photon gauges (I113p . We showed that 
the vacuum polarization has no dependence upon the electromagnetic gauge 
fixing parameter c, but it depends strongly on the two gravitational gauge 
fixing parameters a and h. The effect of being in a general covariant gauge is 
to rescale (155|) by the function K(a,b) given in equation (I124p . By varying 
the constants a and b, one can make K (a, b) assume any values from plus 
infinity to minus infinity. 

Such massive gauge dependence would seem to invalidate any physical 
inference from the results of section 3, however, the gauge independent result 
of Bjerrum-Bohr suggests a simple resolution. There is no question that 
one must include quantum gravitational corrections to the current density 
operator. This seems to be why Bjerrum-Bohr (who implicitly did this) gets 
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a factor of nine different one loop correction to the Coulomb potential. We 
conjecture that making such corrections in a general gauge — which seems 
quite feasible using the techniques of section 4 — would completely cancel the 
gauge dependence of our result. If this could be demonstrated then it would 
be possible to realize the old dream (HI 05] of using solutions to the effective 
field equations as freely as one does classical solutions. Note also that it 
would provide an important class of observables in cosmology, for which the 
S-matrix does not exist. 

The point of this exercise has been to establish the flat space correspon- 
dence limit for a planned investigation of the effects of inflationary gravitons 
on electromagnetism. Our model has been a similar study of the effects 
of inflationary scalars on gravity [23, [21], the flat space limit of which [22] 
played a crucial role in guiding the analysis. In retrospect, we can recognize 
the simplicity of flat space as the ideal venue for sorting out the trouble- 
some issues of dependence upon the choice of field variable and the choice 
of gauge which are so important to a correct interpretation of the many so- 
lutions which now exist to linearized effective field equations on de Sitter 
background [201 EH E21 SSI SZ1 SHI SHI EO] - 
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